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Abstract 

We study entanglement Renyi entropies (EREs) of 1 + 1 dimensional CFTs with classical gravity 
duals. Using the replica trick the EREs can be related to a partition function of n copies of the CFT 
glued together in a particular way along the intervals. In the case of two intervals this procedure 
defines a genus n — 1 surface and our goal is to find smooth three dimensional gravitational solutions 
with this surface living at the boundary. We find two families of handlebody solutions labelled 
by the replica index n. These particular bulk solutions are distinguished by the fact that they 
do not spontaneously break the replica symmetries of the boundary surface. We show that the 
regularized classical action of these solutions is given in terms of a simple numerical prescription. 
If we assume that they give the dominant contribution to the gravity partition function we can 
relate this classical action to the EREs at leading order in Gat. We argue that the prescription can 
be formulated for non-integer n. Upon taking the limit n — > 1 the classical action reproduces the 
predictions of the Ryu-Takayanagi formula for the entanglement entropy. 
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I. INTRODUCTION 

Entanglement entropy (EE) is a powerful observable for many-body quantum systems. 
This is especially so when defined with respect to the reduced density matrix associated to a 
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spatial subregion A of the full system [HE]. EE then detects spatial quantum correlations in 
a fixed many-body state. One simple reason for the appeal of EE is the universal nature of 
its definition allowing for model independent characterizations of many-body phases. To list 
a few applications: EE has been used as an order parameter to distinguish trivially gapped 
phases from those with topological degrees of freedom [3j 0], as a c-function on CFTs in two 
and three dimensions [HH7] and as a measure of thermalization in non-equilibrium situations 
®. 

Unfortunately EEs are rather hard to compute theoretically even for free theories. Tech- 
niques for CFTs are available [9TTTT] and give results for fairly simple spatial regions A. 
However a more general understanding of EE in QFT is lacking. 

Surprisingly there is a simple formula for computing EE in AdS/CFT given by a pre- 
scription of Ryu and Takayanagi (RT) [T2HH] involving the area of minimal surfaces. The 
formula applies to quantum field theories with dual classical Einstein gravity descriptions. 
Some higher derivative corrections have been attempted, see for example [15J, while bulk 
quantum corrections are unknown. The status of the formula remains as a further conjec- 
ture above and beyond the usual rules of the Maldacena conjecture [T6HT8] . In principle one 
should be able to derive it using just these rules, however the attempt in [19] failed as was 
emphasized in [20]. In particular a derivation would forge the way to understanding bulk 
quantum and classical corrections to the formula. 

The focus of this paper will be 1 + 1 CFTs where the sub-region of interest A is the union 
of a set of intervals along the spatial axis [TOj [21] and we consider only the vacuum state of 
the CFT. The RT prediction for this case was discussed in [20], [23] and involves the lengths 
of bulk geodesies which we summarize in Figure [TJ We will attempt to prove the RT formula 
for this case using the replica trick. This trick involves calculating the Entanglement Renyi 
Entropies (ERE) as an intermediate step 

S n = ^lnTr^r (1-1) 

n — 1 

where is the reduced density matrix in the vacuum of the CFT for the Hilbert space 
associated to the intervals A. The EREs are defined for integer n > 2 and can be calculated 
by the partition function of the CFT on a surface Ai of genus n — 1. Assuming one can 
analytically continue the partition function to non-integer n then the limit n — > 1 gives the 
von Neumann entropy expression for the EE. 

This paper was inspired by some of the results of Headrick in [20J where the ERE for two 
intervals and n = 2 was found for CFTs with gravitational duals. We attempt to generalize 
Headrick's results by finding the gravity solutions which are needed to compute the EREs 
holographically for n > 2. We seek handlebody solutions whose conformal boundary is the 
genus n — 1 surface Ai. To generate such solutions we need to represent Ai in terms of 
its so called Schottky uniformization. This representation of A4 can be roughly described 
as a connected domain in the complex plane with certain identifications on the boundaries 
of the domain. Schottky uniformization allows us to find the bulk handlebody solution 
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FIG. 1: The RT prescription for computing the EE in 1+1 dimensional CFTs for 3 disjoint intervals. 
The CFT spatial direction is a and r is the radial direction of the dual AdS^. The minimal surfaces 
are simply geodesies connecting the ends of the intervals. The sum of the regularized lengths of 
these geodesies computes the EE. There is more than one minimal set of such geodesies and one is 
instructed to find the global minimum. We have shown only 2 cases out of a total of 5. 

by extending the domain boundaries and identifications into the bulk radial direction in a 
particular way. 

Actually there is an infinite set of such gravitational solutions. At finite Newton's constant 
Gn one expects all of these to contribute to the partition function as 

Z M = £)exp(-5£ + 0(G%)) (1.2) 

7 

where S^ r oc G^ 1 is the gravitational action for the classical solution labelled by 7. However 
in the classical limit where Gn — > only the least action solution will dominate and we only 
need to find this one. In this paper we show that one can easily construct a small finite 



subset of the infinite set of solutions that contribute to (1.2). 

Interestingly the solutions we can construct in this way have the property that one can 
formulate a simple numerical problem which computes their gravitational action. The answer 
can then be found numerically for integer n > 2. This formulation can be continued in the 
replica index n to non-integers. This is true despite the fact that the bulk solutions no 
longer make any sense. The limit n — > 1 can be studied exactly and the actions computed 
in this way reproduce the Ryu-Takayanagi formula for the EE involving the lengths of bulk 
geodesies. 



Unfortunately since it is more difficult to construct the missing solutions in (1.2) to check 
that they are all subdominant we are left only with a partial result. The gravitational 
actions we compute via the numerical prescription can only be related to the EREs if we 
assume they are in fact the dominant ones. If one could show that this assumption is correct 
then we could compute the EREs and prove the RT formula. A simple way to characterize 
the missing handlebody saddles is by the fact that the bulk solution breaks some of the 
symmetries of the boundary manifold including for example the replica symmetry which 
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interchanges the different replicas. 1 This would be an interesting phenomena if it were to 
happen, and the investigation of this possibility is left to the future. 

Although we will discuss some results for multiple intervals most of the discussion will 
be for the case of 2 intervals. We expect our results to generalize to multiple intervals. 

Our results match the calculations of a complementary paper [22] which takes the CFT 
perspective to this problem. CFTs with large c and a a small number of low dimension pri- 
mary operators were considered. The arguments in [22] are based on semiclassical conformal 
blocks. We comment more on this paper in the discussion. 

The paper is organized as follows. In Section II we introduce the replica trick which 
tells us to compute the partition function of a certain genus n — 1 Riemann surface A4 the 
properties of which we also discuss here. In Section III we give a numerical prescription 
for computing EREs in 1 + 1 CFTs with a classical gravity description. This prescription 
remains a conjecture since we could not rule out the possibility of other saddles being 
dominant. However since the end result we found is rather simple it is useful to present this 
before delving into the details of its derivation. We subsequently show that these saddles 
reproduce the RT prescription and reproduce other known results in the literature. In 
Section IV we discuss the essential ideas behind the program of Schottky uniformization. In 
the Section V we gives details of the bulk solutions that we find. In Section VI we compute 
the bulk action in a few ways and relate the answer to the prescription given earlier on. We 
end with a discussion. There are several appendices with details. 



II. THE REPLICA TRICK AND THE RIEMANN SURFACE 



We are interested in computing the ERE for a spatial region A - the set of N intervals: 

A = [zx, z 2 ) U [z 3 , z 4 ] . . . U [z 2N -i, z 2 n] (2.1) 

where the are cyclicly ordered. The Hilbert space factors locally: H = H A ® H^c where 
A c is the complement region to the above intervals. The EE in the vacuum state is defined 
by: 

p^ = Tr^ c |0)(0| -> S EE = -Tr A (p A \ogp A ) (2.2) 



and the ERE generalizations were given in (1.1) such that lim^i S n = See- The replica 



trick allows one to formulate Tr(p^) n as a partition function of the theory on a particular 
manifold. The arguments are standard and can be found for example in the review [21] . For 
each of the n factors of p A one introduces a Euclidean path integral on the complex z-plane 
with certain boundary conditions on the z real axis. The trace and sum over intermediate 



There are also non-handlebody solutions which are usually assumed to be sub-dominant since they would 
be pathological from an AdS / CFT point of view. We come back to these as well as the replica breaking 
saddles in the discussion. 
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states then glues together these n copies of the z-plane along the intervals in A in a particular 
way. The result is a Euclidean path integral on an n-sheeted Riemann surface or branched 
covering defined by: 

i=l...N v ^ y 

with the entanglement region .4. lying on the real z axis. 2 The genus of this surface is 



(N — l)(n — 1). Beyond this specifying (2.3) only tells us the complex structure of the 



surface, however to compute the CFT partition function we also need to give a particular 
metric in the fixed conformal class. We take this to be the original metric that the CFT 
lives on: 

ds 2 = dzdz (2.4) 

On the branched covering this metric necessarily has conical excess singularities at the branch 
points. These can be resolved by cutting out a region e from the branch points and replacing 
the singular metric with a smooth one. The details of this procedure are standard and given 
in Appendix [Cj The Euclidean path integral on M. can be used to compute entanglement 
Renyi entropies: 

S n = — (\nZ M (ds 2 ) -nlnZi) (2.5) 

71—1 

where Z\ is the partition function of the theory on the flat z plane without any branch 
points. 



The isometries of the surface (2.3) include Z n cyclic rotations of the replicas and the 
anti-holomorphic involution which reflects about the real z axis (the symmetry associated 
to complex conjugation due to the fact that the z% all lie on the real axis.) Together these 
generate the dihedral group D n and we refer to this as the "replica symmetry". For more 
discussions on the relevance of these symmetries to computations of the ERE see |24j . 

It is common to think of Zj^ as the correlation function of twist operators in the product 
orbifold theory of n copies of the CFT under consideration: 

Z M oc (ai(zi)a-i(z 2 ) ■ . . ct^n-^ct-i^n)) (2.6) 

up to some regulator factors which deal with the divergences associated to the conical 
singularities. The twist operator o\ enacts the generator of cyclic permutation of the n 
CFTs upon circling it. And the operator cr_i acts inversely to a\. See for example [37] 
whose results are relevant for computations of EREs for general CFTs. The dimension of 
these twist operators is fixed by the central charge c of the CFT: 



2 Some properties of this surface for two intervals N — 2 are summarized in Appendix A. 
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III. PRESCRIPTION 



In this section we give a prescription for finding and computing certain saddles of 3 
dimensional Einstein gravity that contribute to Zj^ by the usual rules of AdS/CFT. Many 
things will be introduced in an ad-hoc way leaving their justification to later. We also leave 
discussions of the explicit bulk solution to later sections. 

The prescription reproduces several known cases as well as the Ryu-Takayanagi formula. 
Throughout this section we will assume that one of the saddles we construct is the dominant 
solution, and thus at leading order in 1/Gjv computes the ERE. It should be kept in mind 
that this might not be the case. And so the prescription given here remains to be proven. 

We claim that in order to compute S n holographically one should use the following recipe: 

1. Consider the ordinary differential equation (ode) defined on the Euclidean z-plane 
with the points Z{ lying on the real axis: 

^(*) + k„tK*)=0; T zz = (r^-2+^-) (3- 1 ) 

where A = (n 2 — l)/(2n 2 ). The pi are called accessory parameters. 



2. Tune pi such that the solutions of (3.1) have trivial monodromy around a set of N 



cycles Cm in the z-plane with the points Zj removed. We label this set by, 

T = {C M :M = 1,...N} (3.2) 

The Cm are defined to be simple non-intersecting (homo logically) independent and 
non-trivial and each encircle an even number of the z^. At fixed N there is some 
number Mn of independent configurations of cycles T 7 which we label by 

T N = {V 1 : 7 = 1,... AM (3.3) 

3. For a fixed configuration T 7 G T/v the monodromy conditions determine the pj. From 
these construct the following "saddle" Renyi entropies S% by integrating: 

= ^—pl (3.4) 

where c is the central charge. 

4. The true ERE is claimed to satisfy: 

S n = min (3.5) 

7 

We give some clarifying comments: 
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The solution ip(z) will later be used to construct a bulk gravitational solution. 

Prior to imposing the monodromy conditions the accessory parameters are real and 
unconstrained except for the three conditions: 

i i i i 

such that the point z = oo is not a singular point of the ode. Thus the point z = 



oo has trivial monodromy and one can think of (3.6) as being contained within the 
monodromy conditions on the cycles in T 7 . 

The counting of the number of unique configurations of cycles proceeds recursively. 
As we add one more interval iV — 1 — > N we can use configurations Tjv-i to construct 
those in Tn- This is illustrated in Figure [2] 




FIG. 2: A recursive argument for generating configurations of cycles in Tn- The black solid lines are 
new curves. The other curves are represented by the shaded blob and are taken from a configuration 
T G T/v-i or Tn-2 as indicated. The last term subtracts off some over counting of the previous 
two terms. The answer is Mn = 3A//v-i — A/jv-2- 



• The condition that each cycle encircles an even number of points Zj is related to the 
fact that these cycles actually live on the Riemann surface M. and we want them to 
come back to the same sheet. 3 

• We will sometimes refer to a given 7 as a saddle since it will ultimately correspond to 
a particular three dimensional gravitational solution. The monodromy conditions on 
Cm G T 7 will tell us which cycles of the manifold M. are contractable inside the bulk 
three dimensional handlebody solution. 

• For a manifold of genus (N — l)(n— 1) we should pick (N — l)(n — 1) non-intersecting 
cycles (out of 2(iV — l)(n — 1)) to be contractable in order to specify a unique handle- 
body. We will sometimes refer to these as "A-cycles." So far we have specified N — 1 



3 Note the non-crossing condition on the cycles is still appropriate despite the fact that some of the cycles 
actually move into the second sheet. This follows from the comment on cyclic symmetry. 
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of these not including one of the cycles in T 7 which is not independent due to the mon- 



odromy condition at infinity in the z plane (3.6). As we will see the remaining cycles 
are related to these by demanding the bulk solution respects the replica symmetry. 
That is we are also implicitly picking a basis of A-cycles: 

{g m (C M ) : m = 0, . . . n - 1; M = 1, . . . N - 1} (3.7) 

where g enacts the cyclic replica symmetry and moves Cm to the adjacent sheet of 
the branched covering. Note that not all of these cycles are independent because 
X]m=o [9™'(C'm)] — 0. This gives the desired (N — l)(n — 1) counting. 

It is easy to see that the anti-holomorphic involution (symmetry under complex con- 
jugation) is also preserved by this choice of cycles. 

Saddles we are missing include ones where the monodromy condition on the Riemann 
surface Ai do not obey the replica symmetry. These cannot be constructed by the ode 



(3.1) which must be generalized in an appropriate way. 



Note that up to some constants T zz will be the expectation value of the stress tensor 
for the associated saddle. It is then clear that A is related to the dimension of twist 



operators (2.7). Furthermore (3.4) follows from applying the conformal Ward identity 



to T zz and comparing to the conformal transformation of the twist operator correlation 



function (2.6) (albeit on a saddle by saddle basis.) We will derive (3.4) later using the 



bulk action for the constructed solutions. 

The central charge is related as usual to the bulk Newton's constant c = 3/(2(7 jv)- 
The prescription above is for large central charge, otherwise the different bulk solutions 



will all contribute to (1.2) including the ones we have not constructed. 



Each saddle will have a counterpart set of geodesies which we can identify with a 
locally minimal surface of the RT prescription. These geodesies can be constructed by 
noting that the configuration of cycles T 7 partitions the Zi into pairs: 

P 7 = {(z i ,z j ) K ;K = l,...,N} (3.8) 

such that (zi,Zj) G P 7 are either both inside or both outside every cycle Cm G T 7 - 
Joining these pairs by geodesies gives the counterpart RT saddle. The homology 
condition which is part of the RT prescription [14J is satisfied for these geodesies. See 
Figure [3] for an example of this. 



As a zeroth order check we consider N = 1 where we find that the conditions (3.6) are 
sufficient to fix the pi, 

2A . . 

Pi = ~P2 = --, r (3.9) 
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FIG. 3: A picture of the correspondence between boundary cycles for a fixed configuration T 7 E 75 
and the bulk geodesies of the RT formula (green curves hanging down from the boundary). The 



geodesies connect points defined by P 7 in (3.8). Notice that in this picture the cycles in T 7 are 



contractable in the bulk without crossing the geodesies. 



There is only a single configuration and it is clear that the monodromy is trivially fixed to 
zero when passing around the points Z\ and Zi- Integrating this we find the standard CFT 
result [1] 

S n [N = 1] = °- ( 1 + 

where e is a UV cutoff kn=i is unfixed, but scheme dependent. 

In general it is a difficult problem to carry out the above steps for N > 1. Firstly 
there is no analytic way to solve for p] given a monodromy condition 7, so one needs to 



ln((z 2 - zi)/e) + k n=1 



(3.10) 



proceed numerically. Secondly one needs to integrate (3.4) to find S%. We assume that the 
partial derivatives commute so we only have a single integration constant for each saddle. 
Clearly the issue here is that the relative integration constants for the saddle entropies are 
not fixed by the above prescription. We will fix these constants by taking limits where the 
saddle entropies S% are related to the results for one less interval. This allows us to fix the 
integration constant recursively. 

Note also that we can give an absolute formula for computing S% which will fix this 



integration constant, see (6.40) for the two interval case. The formula is written in terms 



of the solution ip to the ode problem but is more complicated than the prescription given 
above so we leave that till a later section. The prescription we have given is sufficient for 
our current purposes. 



A. Reproducing the Ryu-Takayanagi prediction 

We wish to compute the monodromy matrices in the replica limit 5 n = n — 1 — » 0. We 
can do this using perturbation theory. Assume that p] vanishes linearly in the replica limit: 
p] ~ pi5 n + 0(5%) where pi are constants which we need to determine. The fact that p] 
should vanish in the replica limit follows from the conjectured formula for the entanglement 



Renyi entropy (3.4) which we expect to be finite in this limit. Also note that A = 5 n + 0(5^) 



in this limit. The second order ode (3.1) can be conveniently represented as a first order 



10 



system: 



The monodromies are then simply path ordered exponentials: 

u{z) = M{C)u{z ) M{C) = Pexp ( / dzfffc) ) (3.12) 



\JC 

where C is a specific path from ;?o to z. Note that M(C) has unit determinant which follows 
from the Wronskian condition of two solutions to the ode. Perturbatively we have 

So we can then use time dependent perturbation theory methods to solve this problem where 
z is thought of as "time" . Firstly move to the interaction picture: 

M(C) = M (z)M I (z) ^- z M I (z) = 5 n (M^ 1 H l M )(z)M I (z) (3.14) 

where the zeroth order solution to the ode ipo = Az + B can be used to find the zeroth order 
monodromy matrix: 

M °M=(d {z T ] ) (3 ' 15) 

The path ordered exponential expression for Mj can then be computed to first order by 
simply expanding the exponential: 

Mj « 1 + 5 n [ dzM^HM = 1 + 5 i [ dz ( {Z ~ Zo) {Z ( ~ Zo)2 ) Tt (3.16) 
Jc 2 J c y 1 -{z-z )J 

If we close the cycle C by sending z — > Zq we find the monodromy condition requires the 



vanishing of resulting contour integral in Eq. 3.16 For a cycle C = Cm £ T 7 we get three 
independent conditions: 

^p, = 0, ^(p^ + l)=0, ^(p^ 2 + 2^) = 0, (3.17) 



where the sum is over points Zi contained in the interior of Cm which we have denoted by 
the domain Dm- Note that it does not matter which "interior" we choose - because of the 



monodromy condition at oo given in (3.6). After some thought it becomes clear that this 
set of N equations is solved by the following conditions on the pairs (zi, Zj) G P 1 into which 
the cycles Cm partitioned the Zj. 

Pi = — , Pj = ^— V Zj ) e P 1 (3.18) 

Z% Zj Zj Z{ 
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The saddle entanglement entropy can then be found by taking lim n _ s> i in (3.4) and integrating 
the result: 

<^ = f E H\^-Zi\/e) + Kl (3.19) 

This result is exactly c/6 times the regulated lengths of geodesies in AdS^, connecting the 
points (z;, Zj) e P 7 on the boundary. As we discussed around Figure [3] there is a corre- 
spondence between the saddles that we construct (at any integer n > 2) and the minimal 
surfaces (geodesies) needed to compute the RT answer. We have shown here that the action 
of these saddles can be continued to non-integer n and in the limit n — > 1 they become the 
lengths of the corresponding RT geodesies. While other aspects of this section are somewhat 
conjectural, the last statement is correct and hints at the inner workings of the RT formula. 

To completely reproduce the RT prescription we are left to compute the integration 
constants k 7 relative to all the different saddles T 7 G T/v- We give the following argument. 
Firstly consider an adjacent pair (zk,Zk+i) € P 7 which is enclosed by a unique single cycle 
Cl € T 7 which does not enclose any other z% . Note that there must be at least one such 
pair. Now take the limit Zk —> Zk+i where we expect to reproduce the entanglement entropy 
for N — 1 intervals and a configuration of cycles given by (r 7 / = F 7 \Cl) G 7jv-i- Atleast 
up to a UV divergence associated with the closing of the interval [zk, For a very small 

interval Zk ~ Zjt+i we can zoom in on this and ignore all the other intervals - allowing us to 
exactly subtract off the EE associated with this single interval. Note that it might be that 
[zk-, 2fc+i] is not an interval in A but is an interval in the complement A c . In which case we 
can appeal to approximate purity of the state at small distances so this still contributes the 
same divergence. Further we require that we are in a regime where S EE is the dominant 
saddle - this should be possible to arrange for by moving around the other Zj. We find that 

lim [S EE - I ln((zjb+i - z k )/e) - k ) = S]^ E (3.20) 

Note we are assuming that the regulator we use is such that it treats the UV divergences 
located at the different points Zi in a uniform way. This way we get exactly S EE on the right 



hand side of (3.20) and no other ambiguous constants. Assuming that «y = (N — 1)k\ is 
fixed for all configurations in 7jv-i then by induction we find iC N = Nk,i which also must 
hold for all saddle configurations in Tn- The final answer: min 7 S EE is then the RT formula 
for disjoint intervals in a 2d CFT. 



B. Two intervals and the mutual information 

We now specialize to the case of two intervals N = 2. We think that most of the following 
results work for N > 2 but the arguments become cumbersome and we content ourselves to 
looking in more detail at the first non-trivial case. According to our prescription we have 
two different configurations of cycles which we label 7 = a, /3. These cycles are shown in 
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Figure |4} For example they correspond to the following partitioning of the z% into pairs: 

P a = {(z 1 ,z 2 ),(z 3 ,z 4 )} , Pp = {(z u z 4 ),(z 2 ,z 3 )} (3.21) 

For ease of notation we will often drop the 7 = a, /3 subscript when the distinction is not 
important. 






FIG. 4: The case for two intervals. The black sold curves show the monodromy cycles V £ T2. 
There are two sets of cycles which we label a,/3. The dashed curve on the right figure is due to 
the fact that this curve actually moves through the branch cut into the last (n'th) replica. We also 
define here the specific cycles C a , Cp for later reference. 

To start we would like to understand more about the dependence of the prescription on 
the Z{. The conformal transformations that leave the z plane vacuum invariant and also 
leave the Zi on the real axis are given by SX(2,R) transformations: 



Az + B 
Cz + D 



Zi -> 



Az, L + B 
Czi + D 



AD-BC = 1 



(3.22) 



where A,B,C,D are all real. These move around the points z-i and do not change the 
ordering up to cyclic permutations. We would now like to track the transformation property 



of Pi in (3.1) under SX(2,R). Firstly the conditions (3.6) which the Pi satisfy are left 



invariant if we transform: 



Pi ->■ {Cz % + Df (pi + 2A 



C 



Czi + D 



(3.23) 



This property makes it clear that pi transforms almost like a differential d z .. More precisely 



according to (3.4) pi is conjugate to Zi and so in order to reproduce the transformation 
given in (3.23) we must demand that the entanglement entropy is not SL(2, R) invariant 
but rather: 



S n $n 



6 V 1 + n) ^ 

i=l 



4 



; ln(C^ + D) (3.24) 
These statements follow trivially from the fact that ERE can be represented as the log of a 



four point function of twist operators (2.6). Then the SX(2,M) transformations above are 
simply due to the conformal weights of the twist operators. 
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In order to soak up this transformation we can define what is known as the Mutual Renyi 
Information (MRI). The MRI is the following combination of entanglement Renyi entropies: 

I n = I n ([z 1 ,z 2 ), [z 3 ,z 4 ]) = S n ([zi,z 2 )) + S n ([z 3 ,z 4 )) - S n ([z 1 ,z 2 ] U [23,24]) (3.25) 

such that: 

In = -S n +^(l+ ^ In ((z 2 - Zl )(z 4 - z 3 )/e 2 ) + 2k, (3.26) 
It follows that I n is SX(2,R) invariant and as such can only depend on the cross ratio: 

{z 4 - z 2 )(z 3 - Zi) 

From this we can define the following SL(2,M) invariant accessory parameter: 

dx ~ 6(n-l) Px [6 - Z * } 

There is one final property we have not exploited. Since we are working in vacuum (to 
define our density matrix pX) it follows that the Renyi entropies satisfy = S^c where A c 
is the complement of the region A. This implies the following: 

S n ([ Zl , z 2 ) U [z 3 , z 4 )) = S n ([z 4 , Zl ) U [z 2 , z 3 )) (3.29) 

This purity relation corresponds to switching z 4 -h- z 2 and as such can be thought of as a 
very simple crossing relation for the twist operators. When we plug this into the mutual 
information we find: 

I n {x) = 4(1 - x) + l(l + ±) In (j^A (3-30) 

Actually we can go a little further. If we track what happens to the configurations of 
cycles in % as we send z 4 -H- z 2 it is clear that r Q -B- Yp. So in terms of the saddle ERE we 
can define the saddle MRI which satisfies: 

^HW-^ + ^l + i)!*^) (3.31) 

and similarly for a -B- 0. 

In all we can now refine our prescription a little more. We can make a conformal trans- 
formation to move the points to Z \ = 0, z 2 = x, z 3 = 1, z 4 = 00. After which our ode looks 
like: 

, „ , , 1 / A A A 2A p x x(x-l) \ . 



and the prescription to compute the mutual information is simply to integrate (3.28). We 
should then find the maximum of the two possible saddle mutual informations (note the sign 
switch in the definition (3.26) of I n ): 

I n (x)=max{IZ(x),lZ(x)} (3.33) 
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where we remind the reader that it is possible there are some missing saddles which become 
dominant at some x and thus override this answer. 

It is clear that in the limit where x — > the dominant configuration is T a and in this 
limit the mutual information vanishes since this limit corresponds to moving the two intervals 



infinitely far apart. This condition will be used to fix the integration constant in (3.28). 

Since when x = 1/2 the two different saddle mutual informations agree (by purity of the 
vacuum state) it must be the case that the T a and Tp saddles switch dominances at x — 1/2 
[20] . This results in a first order phase transition for any n. 

If we are feeling lazy we can reconstruct the calculated contribution to the mutual infor- 
mation from a single saddle: 

I%{x) , < x < 1/2 

*-M-< W -, ) + |( 1 + IW^_). vx.o (3 ' 34> 



C. Reproducing the known answer for n = 2 

Set n = 2,A = 3/8 and it turns out in this limit we can analytically solve the ode. The 
reason lies in the fact that we are in this case secretly describing a genus one torus. The 
case for n = 2 was already worked out in |20j based on fairly extensive computations given 
in We will see that our prescription reproduces their results with relative ease. 

The two indepdenent solutions can be written as: 

= JiAm ex P (± ht ( z )) > * (*) = 7 , * ( 3 - 35 ) 

where h is an unfixed constant which is related to the accessory parameter: 

2-x h 2 
Px ~ Ax{x - 1) + 2x{x - 1) (3 ' 36) 

These solutions can then be used to find the monodromy matrix: 

( hf c t'(z)dz Q \ 

e-W*) ^ Zo) (3 ' 37) 

where C is a path from zq to z and the non-path dependent factors are: 

The trivial monodromy condition for the curve C G T a is then simply: 

2nik = h f t\z)dz = 2h f dZ = = AhK(x) (3.39) 

Jc Jo v z(z — l)(z — x) 
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where k is an integer and K(x) is the complete elliptic integral (defined by the integral 
above.) The integer k is unfixed so far. The case k = does not work since one has to 



be carefully when taking h due to the degeneration of the matrix ^f(z) given in (3.38) 



in this limit. For \k\ > 2 the solutions one finds involve a multiply wound uniformization 



coordinates (see Section IV) and should not be included in our prescription since they will 
not correspond to sensible bulk solutions. We are left with k = ±1 of which either gives the 
same answer. The accessory parameter is: 



X 7T 2 



Which integrates to: 



Px 4x{l-x) 8x{l-x)K{x) 2 ^ 3 ' 40 ^ 

/« = -^log(2 8 (l-x)/x 2 ) + ^r 2 (3.41) 

where we have defined the (purely imaginary) modular parameter for the underlying torus: 

K(l-x) 

n - = "Wr (3 - 42) 

and we have added an integration constant such that If (x = 0) = 0. Similarly we can find 

1% by imposing the different monodromy condition on the cycles Tp. The answer one finds 



satisfies the expected purity relation (3.31 ) where under x — > 1 — x the modular parameter 
of the torus undergoes the modular S transformation r 2 — > 1/V 2 . This makes sense because 
S switches the cycles on the torus and thus the two monodromy conditions we are working 
with. 



D. Numerics for n > 2 

In order to compute the monodromy matrices numerically it is convenient to define con- 
nection matrices along the real line between the singular points. These matrices relate 
canonically chosen linearly independent solutions at adjacent singular point. We relegate 
the details to Appendix [Bj The monodromy condition can easily be read off from these 
connection matrices and from this we can compute p x . 

The results are shown in Figure [5] Actually there is very little difference between the 
Mutual Renyi Information for different values of n and in order to effectively compare them 
we subtract off a scaled version of the EE (n = 1) which takes into account the scaling of 
the twist operators with n: 

J n (x) = I n (x) + I x [x) (3.43) 

Recall that the RT formula for the MI is I\{x) = max{0, (c/3) In The function J n (x) 

has the property that it is symmetric about x = 1/2. 
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FIG. 5: Calculated contributions to the Mutual Renyi Information (MRI) in holographic CFTs. 



We show in the left panel a subtracted version of the MRI as defined in (3.43). Both saddles T 



and are important and dominate for x < 1/2 and x > 1/2 respectively. In the right we show the 
dependence of the MRI on re for fixed x (after analytically continuing from integer re). The re = 1 
limit for x < 1/2 is zero as predicted by the RT formula. Numerically it was more convenient to 



use (6.40) to find this right plot. 



IV. SCHOTTKY UNIFORMIZATION 



Having introduced the ode (3.1) as the main crux for constructing certian bulk solutions 



we should now explain where this came from, and in particular give some pictures of what 
the bulk solution looks like. Since all solutions of Einsteins equations with a negative cosmo- 
logical constant in 3 dimensions are globally quotients of AdSs we simply need to determine 
the appropriate quotient. The technology we need in order to do this goes under the name 
of Schottky uniformization. We give here a rough general discussion of this technology. We 
follow closely the discussion in [27] and [2B] see also 
Pick the following coordinates on AdS^\ 



ds 



d^ 2 + dwdw 

= e 



(4.1) 



with conformal boundary at £ — > the complex w plane. The isometry group of AdS% is 
PSL(2, C) where the action induces a conformal isometry on the boundary. The action on 
AdS 3 is: 



w —> 



(aw + b) (cw + d) + ac£ 2 



\cw + d\ 2 + |c| 2 £ 2 
where ad — be — 1. As £ — > this action becomes: 



\cw + d\ 2 + |c| 2 £ 2 



(4.2) 



aw + b . . .... . 

ir -»> ; = L(w) £ -»■ £\L'(w) 



cw + d 



(4.3) 
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In this way the quotient of AdS^ by a discrete subgroup £ C PSL(2, C) descends to a 
quotient of the complex w plane. 4 This quotient is then a way of representating the surface 
A4 = C'/E. Thus one of the steps we will need to understand is how to map the w complex 
plane into Ai: 

n s : C ->• M (4.4) 



consistent with the action of the quotient. In fact the ode (3.1) is exactly what determines 



this map. At the same time the monodromies of the ode (3.1 ) determine the correct quotients 
of the w-plane. 

Let's see roughly how this works. Consider instead the inverse map w = n^^z) which 
is multivalued on Ai. The Schwarizian derivative of this map behaves like a holomorphic 
CFT stress tensor on Ai: 

This equation can be thought of as a differential equation for -k^^z) = w(z) where T zz is 
taken as a fixed input. We will construct T zz independently in a moment. Solving equation 



(4.5) is equivalent to solving a second order ode: 



w(z) = ^\, r+ \T„1> = (4.6) 
where ip = ipi2 are two linearly independent solutions of this ode. At this point we have made 



a connection with the prescription of Section III However we still need to give an argument 



that the tress tensor T zz takes the form quoted in (3.1). Since the map n^ 1 is globally 
defined (but multivalued) and since the the Schwarzian derivative does not change under 
the PSL{2, C) action on w the stress tensor is globally defined on Ai and not multi- valued. 
However T zz does not transform homogeneously under general conformal transformations 
z z(z) since the Schwarzian derivative shifts under such coordinate changes: {t, lz} = 
z'(z) 2 {t, z} + {z, z}. The rule on the overlapping patches is: 
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T~ z =[ — \ T zz + {z,I} (4.8) 



4 We are being heuristic here - for example we should first remove a certain set of measure zero from the w 
plane, for which T acts badly (fixed points of T): C = C/{bad points}. We can then form the quotient 
C'/S. For a proper discussion see. We will continue to be heuristic without making similar admissions. 

5 Note that this also implies that the solutions of the ode transform as —1/2 differentials: 



</;(2) = ^(*(2)) ( m ) (4.7) 
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Given these properties we claim that the following expression (on the z coordinate patch) 
is smooth on A4 and completely general: 




-Z 3 + Zi+ z 2 + z 4 
(z - zi)(z - z 2 ){z - 




(4.9) 



where A = 1/2(1 — 1/n 2 ). To argue for this form consider to begin with the the last sum 
over oj s zz which are holomorphic quadratic differentials on Ai. There are 3(n — 2) of these 
and they are enumerated in Appendix |A| Given a fixed T zz we can always add a linear 
combination of quadratic differential since they transform homogeneously under conformal 



transformations thus preserving (4.8) 



Finally we have to check that the remaining term multiplying A in (4.9) is smooth on 
M.. This term is not a quadratic differential which are smooth by definition. We only need 
to check the behavior near z — V Zi and z — > oo. Close to for example z — > z\ we can use the 



coordinate y defined by the branched covering (Al) y n ~ (z — z\) such that the Schwarzian 
derivative is: 

1 n 2 — 1 

{z,y}^-~ — + 0(y n - 2 ) (4.10) 

* y 

and thus the new stress tensor in the y coordinate patch is: 

T yy (y) = j % (An 2 - \{n 2 - 1)) + 0(y n - 2 ) (4.11) 
This is smooth provided A = 1/2(1 — 1/n 2 ) and n > 2. The last term in the brackets of 



(4.9) is then required for smoothness as z — > oo. 

The logic of the preceding discussions is that we have replaced the problem of finding 
the map 71-^ 1 with the problem of finding the accessory parameters p s . As was discussed 



extensively in Section III these should be determined by the monodromy conditions imposed 
on solutions to the ode. For example if we traverse a closed path ConM the map w = n^ 1 



defined by (4.6) is not single valued since the solutions (t/>i, ip 2 ) undergo a monodromy M(C): 



{^M^{^MM{C) =► uj^^^ ( aC )=M(C) (4.12) 

cu + a \b a J 

Thus the monodromies of the ode determine a PSL(2, C) action on w. And in this 
way they determine the discrete quotient group S. Note that the monodromies form a 
representation of the fundamental group of the surface M. and so does the quotient group 
S. 

We need to understand which groups E produce the desired handlebody solutions when 
acting on AdS^. These groups are called Schottky groups and we simply quote some results. 
They have the property that £ is freely generated by half of the 2(n — 1) generators in 
the fundamental group. We define these generators through their PSL(2, C) representative: 
{L m : m — 1 . . .n — 1}. Upon traversing around the other half of the generators of the 
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fundamental group one finds trivial monodromy and trivial action in the quotient group. 
These generators correspond to a basis of non-intersection cycles in the homology of Ai, a 
basis of "A cycles" . It turns out that these are the cycles which are contractable in the bulk 
of the corresponding AdS% handlebody. The "B cycles" then correspond to the generators 
which have nontrivial action L m and are not contractable. 

We can get a rough picture of the quotient looking more carefully at the fundamental 
domain of C'/S and AdS^/Y^. See Fig. [6] for a picture. The fundamental domain is given by 
specifying (n — 1) pairs of non-intersecting circles in the w plane: {C m , C m : m = 1 . . . n— 1} 
and then identifying the circles C m and C m via the non-trivial generators L m {C m ) = C m . 
Note that the fundamental domain is not unique. The generators L m map the outside of the 
circle C m into the inside of the circle C m . So for example we can shrink C m while making 
C m larger and still have a fundamental domain for the quotient. 































Q, 















FIG. 6: The Schottky fundamental domain for a genus 2 surface with 2 generators L\,L2- The 
shaded region is the domain (continued to infinity.) Note that typically one normalizes the gener- 
ators using the freedom to conjugate by PSL(2,C) such that one of the circles surrounds w = oo 
which is then absent from the domain. In the above picture we have not done this, since this will 
be convenient for us later. 



The fundamental domain of the quotient of AdS% is simply found by extending the cir- 
cles C m , C m living on the boundary to hemispheres in the bulk of (4.1). These are two 
dimensional minimal surfaces in AdS^ ending on the circles. Note that on M. the cycles 
which encircle C m , C m are contractable within the three dimensional bulk solution. These 
correspond to the cycles of Ai with trivial monodromies (the A-cycles.) The B-cycles are 
paths in the fundamental domain which connect the identified circles. 

Finally we come to the accessory parameters. In (4.9) we have 3(n — 2) of these, however 
we claimed in Section III that there was only a single independent accessory parameter for 
two intervals p x . Note that of all 3(n — 2) quadratic differentials enumerated in Appendix |A| 
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only one of them lo\ z does not change under the actions of the replica symmetry. This 
means that in order to preserve this symmetry p s = for s ^ 1. If we include the anti- 
holomorphic involution (complex conjugation on the z-plane) in the replica symmetry we 
find that the remaining accessory parameter pi = p x should be real. Fixing most of the 
accessory parameters to zero is only possible if the monodromy conditions respect the replica 
symmetry otherwise these should be turned on and the resulting bulk solution will also not 
by symmetric. 



V. BULK SOLUTION 

We turn now to a detailed description of the bulk solution, from which the final goal is 
to compute the bulk action which we get to in the next section. We start with the details 
of the quotient C/S. 

Assume the Schottky monodromy problem has been solved. As discussed in the previous 
sections for two intervals and for bulk solutions which are replica symmetric there are two 
different monodromy conditions that we can impose that we labelled T a ,Tp 6 T2 (see Fig- 
ure E]). For arguments sake pick Tp which involves imposing trivial monodromy around a 
cycle Cp which encircles the pair (z 2 , z$). The other case can be worked out in an analogous 
manner. We would like to work out the identification circles C m , C m for the Schottky fun- 
damental domain as well as the generators L m linking them. The Schottky group £ is only 
defined up to common conjugation by PSL(2, C) and thus we can choose two independent 
solutions of the ode at will in order to produce the map w(z). We pick the solutions, 

V>± = (z - ^i) 1/2±1/2/n (l + 0(z - Zl )) (5.1) 

which diagonalize the monodromy around the point Z\. Then define: 

w = A — (5.2) 

V- 

such that w(z\) = 0. Note that under this choice the Z n replica symmetry is generated by 
rotations of the w plane by an angle 2ir/n. A nice way to get a concrete picture of the map 



generated by (5.2) is to consider the images under w(z) of the 4 real axis segments in the 
z-plane between the points Zj. Segments slightly above and slightly below the real z-axis 
±irj should both be considered since these will map to different curves in the w plane. We 
should also consider the segments on all the n replicas. The images of these segments will 
then trace out a particular fundamental domain in w. The identifications L m can be worked 
out by appropriately glueing the real line segments together amongst the different replicas. 
See Figure [7] for the resulting picture. This can be confirmed numerically by plotting 



(5.2) after one has imposed the monodromy condition. We will not give the full detailed 
argument that leads to this picture. However we summarize some of the more important 
aspects: 
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LHP 



x s 1 



FIG. 7: A picture of the map from the branched covering z (left) to the Schottky domain w (right) 
for n = 3. We only show the image of the first replica. The solid line tracks the real line segments of 
the z plane which are slightly above the real axis and the dashed line tracks the segments just below 
the real axis. The blue fuzzy lines represent the branch cuts. The point z\ maps to w(zi) = oo. 
Note that we have imposed the monodromy condition on the cycles in Tp (Figure |4]) so that the 
map w(z) jumps discontinuously across the segment [£2,23] via the generator L\. 



Since the ode (3.1 ) is real along the real z-axis all the segment images must be circular 
arcs or straight lines. That is there is always a basis of solutions to the ode ipi,ipu 
which is real along a given segment. Then: 

a£ + b „ 4>i(z) m 
w seg = —. — ; £ = e R (5.3) 

g c£ + d if> n (z) v ; 

for some a, b,c,d G C. This describes a circle or line in the complex w-plane with 
affine parameter £. 

Most of the real axis segments map to straight lines along rays emanating from the 
origin in the w-plane w = £exp(i2mTr/n) and w = £exp(i(2m + l)ir/n) where m = 



l,...n. In particular w(z±) = 00. This behavior for (5.2) only follows once the 
monodromy condition is imposed. 

Images of the segments z G \zi , 23] ± if] on all the different replicas map to circular arcs 
which meet the above straight lines at an angle of rr/n. These arcs can be described 
as: 

— e -j^~ J 0<£<oo z 2 <a < z 3 (5.4) 

1 — £e n J 

/ \ I xq — £e^ \ , 
w(o - it]) = e tzn n ^— < £ < 00 z 2 < < z 3 (5.5) 

\ I — £e^ J 

where m labels the images generated from the different replicas m — 1, . . . n. 
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We have chosen the magnitude of A in (5.2) so that |iu(£3)| = 1. Then |tu(^2)| = 



is the single remaining parameter which can be computed numerical in terms of the 
cross ratio x. It satisfies xs < 1. 



Note that (5.2) jumps discontinuously across [z%, 2&] on the z-plane since the arcs (5.4) 



are different from (5.5). This jump is encoding a non-trivial monodromy element and 
occurs on one of the B-cycles. Note in particular this jump does not occur across a 
branch cut on the z-plane. 

We can compute the associated generator in S by finding the PSL(2, C) transforma- 



tion which identifies an arc in (5.4) with the corresponding arc in (5.5) (with the same 



m.) By symmetry the two affine parameters i map onto to each other. One finds: 

(■ 27T • (2m — 1) N 

x s -e 1 ^ 2ix s e ln ^r- sm(ir/n) . 

• (2m- 1) . . . -2tt l^-^J 

-2ie~ m -^ sin(vr/n) x s - e~ 1 ^ 1 

• Gluing the replicas together we get the global picture on the left side of Figure [SJ 
Note that the circular arcs are the boundaries of the fundamental domain and they 
are identified pairwise as in the Figure. 

Note that the arcs start to cross unless we demand cos 2 (7r/n) < xs- This actually 
corresponds to the boundary of moduli space x — 0. See FigireflOlfor numerically calculated 



plots of xs(x). The generators (5.6) satisfy 



TrL m = 2 ( ^- C ° s(27r/ " )) (5.7) 
1 - x s 

and the condition that the arcs do not cross also corresponds to the requirement that the 
elements are loxodromix (| TrL| > 2.) If we parameterize L m as: 



L m (w) —a m fw — a 



L m (w)-r m \w — r r , 

then 



(5.8) 



x s - cos(27r/n) - sin(27r /n)y/(x s / cos 2 (vr/n)) - 1 

q = = (5.9) 

Xs — cos(27r/n) + sin(27r/n)y (xs/ cos 2 (n/n)) — 1 

Note the parameter q does not depend on m. The attractive and repulsive fixed points 

r,w. arg(a m ,r m ) = 7T -^- — ± tan -1 (^J{x s / cos 2 (vr/n)) - (5.10) 

The two fixed points come together at the boundary of moduli space when cos 2 (7r/n) = xs- 
The generators are not all independent: L n . . . L 2 Li = (— l) n where —1 acts trivially as 
a fractional linear transformation. This leaves n — 1 generators for the group S. Which is 
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FIG. 8: Pictures of the fundamental domain of the Schottky quotient. We have drawn the case 
n = 3 for some fixed xs{x). The left plot shows a non-standard domain D s which is however clearly 
Z n symmetric. There are three (generally n) generators Li,L2,L% for this case which are shown 
as green arrows. These generators are not all independent. In the right picture we have deformed 
some of the circles in the left picture to form a more standard fundamental domain T>d where there 
are now only two generators L\,L2- Note that the images of the points w(zi) only appear once 
within this fundamental domain. In the right figure we label the blue circles 0\,0\ and the red 
ones 02,02 • 




FIG. 9: The symmetric fundamental domain D s for n = 6. 
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FIG. 10: The single parameter xs which goes into the Schottky group for the surface M. with 
monodromy conditions Tg. For x — > (where this saddle is subdominant to T a ) the limit is 
xs = cos 2 (ir jn) - although seeing this numerically requires high precision. 

the number that we expected. The fundamental domain is pictured in the left of Figure [8] 
for n = 3 and Figure [9] for n = 6 and uses the real z-axis segments that we found above. 
This domain also does not conform to the standards of the usual fundamental domain of 
a Schottky group. This fact goes hand in hand with the over counting of generators. It is 
easy to see how to fix this. By deforming the circular arcs in an appropriate way the last 
generator L n becomes superfluous and one is left with 2(n — 1) identified closed circles rather 
than arcs. The argument is sketched in the right of Figure [8j The existence of this deformed 
fundamental domain for the group generated by L\, . . . L n _i means that this group is by 
definition a classical Schottky group [26]. In what follows we will go back and forth from 
considering these two different fundamental domains. We will refer to the first symmetric 
domain as V s and the later deformed domain as T)^. 

The boundary of the fundamental domain is defined as: 

n n—l 

&D S = |J (U m U U n }j 8V d = |J (0 m U O m ) (5.11) 

m=l m=l 



where U m ,U m are the arcs given in (5.4), (5.5) respectively. While the circles O m ,O m are 



not uniquely defined; an example is given in the right side of Figure [8j 

The attractive and repulsive fixed points all lie along the same circle \w\ = y/x~s in the 
complex w plane. This demonstrates a fact that was speculated upon in [21]. The authors 
showed that if the Schottky parameters a m ,r m ,q m can be chosen to be real then there is 
a so called real duality between the compact boson CFT at the self dual radius and free 
fermions (with a fixed spin structure) and this implied the EREs of these two theories were 
the same for two intervals. They demonstrated the real duality using a different method and 
speculated that this meant one could choose the Schottky parameters to be real. Indeed we 
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see here explicitly that this is the case - by making an PSL(2, C) transformation to send 
the circle \w\ = y/xs to the real axis then a m ,r m will all become real, q remains fixed under 
this transformation and is real as above. 6 

Note that this reality argument applies to the Schottky parameters for our specific choice 
of A-cycles (cycles with trivial monodromy). It probably does not apply for other replica 
symmetry breaking choice of A-cycles, although for the arguments in [24J one only needs 
the reality condition for one such choice of cycles. 

Finally a picture of the bulk solution can be drawn by extending the circles into hemi- 



spheres in AdSz (4.1). We depict in Figure 11 the symmetric case where the fundamen- 



tal domain on the boundary is T> s . In this picture two bulk hemispheres intersect over a 
geodesic. We speculate that one can identify this with a generalized version of geodesic in 
the RT prescription. 




FIG. 11: The gravity solution found by extending the arcs U m , U m which live on the boundary of 
AdSs into hemispheres inside the bulk and identifying these hemispheres. We show again the case 
n = 3. The picture also represents a (non standard) fundamental domain for the quotient of AdS% 



by the action of S which acts as (4.2). The hemispheres intersect over the green curves which are 
bulk geodesies. 



6 This also means that the Schottky group is actually a Fuchsian group acting nicely on the disk \w\ < ^fxs- 
Interestingly this allows us to find a real-time three dimensional black hole based upon M. , see [3§1 130] for 
details. These are generalizations of the usual BTZ black holes [31] ■ It would be interesting to understand 
what such a solution means for the EREs. 
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VI. BULK ACTION 



We will consider two ways to calculate the regularized on shell Einstein action for the 
bulk solution. Certain results in the literature will be used heavily. To begin with we work 
on the deformed domain T>d defined in the previous subsection. It was shown in [27] and 
further in [32] that the on-shell action can be written in terms of a certain two dimensional 
Liouville action living on the boundary. The action is defined on the domain T>d- It was first 
written down in [28] by Zograf and Takhtajan (ZT) and we will refer to it as the ZT action. 
It was further shown in [28] that the variation of the ZT action with respect to the moduli 
of A4 behaves nicely and we will use these results to prove the assertion of the prescription 
that the variation of the EREs gives the accessory parameters (3.4). 

After this we will go through a re-derivation of the results in [27] for the bulk action using 
the symmetric fundamental domain T> s which is somewhat more convenient for our purposes. 
This will allow us to give an absolute expression for the EREs not involving derivatives with 
respect to Z{. 



A. Zograf- Takhtajan Action 

Firstly we introduce the notion of Fuchsian uniformization. We only need it as an inter- 
mediate step so we will be brief. This is another kind of uniformization compared to the 
Schottky variety which aims to place a constant negative (for genus (n — 1) > 1) curvature 
metric on M.. The method is very similar to the Schottky case. Consider the Poincare disc 
D with metric: 

_ . o dtdt j _ , 

where \t\ < 1. Fuchsian uniformization represents M. as a quotient of D by a discrete group 
which acts nicely on it. That is a discrete subgroup of SX(2,C) which leaves the metric 



(6.1 ) invariant. 

Several results in the literature are available for computing the gravity partition function 
when the metric is taken to be ds 2 . Because of the conformal anomaly the result does depend 
on which metric we use within a fixed conformal class. We actually want the partition 



function on ds 2 given in (2.4) and the difference between these two is given by the Liouville 
action: 

Z M [ds 2 } = e s ^Z M [ds 2 ] , S L = ^ J dH (W) 2 - {i m tL 2 ) (6-2) 

where <$f is the conformal factor which relates the two metrics: 

ds 2 = e'^ds 2 = e~ Mz) dzdz (6.3) 

The metric ds 2 has conical singularities which means 4>f is singular at these points. We 
deal with this by cutting out holes around these points which is a standard [47] procedure. 
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Details are given in Appendix [Cj To find (pp we need to map the branched covering to this 



representation of Ai. Once again the ode (3.1) allows us to construct the analytic map and 



the quotient group Si? as the monodromy group: 



where the stress tensor Tf z has the same form as the Schottky case (3.1) however now the 
accessory parameters pi — > pf will be different. In order to fix the pf we must impose 
the condition that all monodromy elements generated by the fundamental group leave the 



metric (6.1) invariant. This is called the Fuchsian monodromy condition. Note there is 



a unique condition here, we do not need to pick different "A cycles" and "B cycles". We 



emphasize that this is a different monodromy problem to the one given in Section |III| and 
that we expect to find a different accessory parameter. The field (j)p is then: 

<j) F = -ln|t'(z)| 2 + 21n(l- \t(z)\ 2 ) (6.5) 

and we are now in a position to compute Sl- 

To calculate Zj^(ds 2 ) we introduce the ZT action which also happens to be a Liouville 
type action however now living on the Schottky w space. Firstly one introduces a new 
Liouville field <$$ which is the conformal factor on the Schottky w-space which uniformizes 
that space, placing on it a constant negative curvature metric: 

ds 2 = e'^dwdw (6.6) 
<Ps = -ln|t'(w)| 2 + 21n(l- \t{w)\ 2 ) (6.7) 

A major difference between the two Liouville fields is that (f)$ is not single valued on 
Ai where as <pp is. As we move around on Ai the w coordinate undergoes PSL(2, C) 
transformations along the B-cycles: w = L m (w). In order that <ps is consistent with these 
jumps in w (the action of S) we must also require that (fis jumps: 

<Ps(w) = <PsH+log\L' m (w)\ 2 (6.8) 



under which the metric (6.6) is preserved: 

e'+sWdwdw = e-tsMdwdw (6.9) 

The field (fis is uniquely specified by the identifications (or boundary conditions) and the 
requirement that it satisfies the Liouville equation: 

dA<f> 8 = -2e-* s (6.10) 



This equation follows from the requirement that (6.6) has constant negative curvature. 



More succinctly <$$ is the solution of the equations of motion which follow from varying the 
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following ZT action denned in 



96vr 



/ d 2 w ((d<f) S ) 2 + 16e-* s ) + S 

J-Dd 

lb J. n 

SzT = J2 / (205 + log|^| 2 + 21og 

m=l 



96tt 



ra-1 



r II tii 

I2A / j^^i _ idw^ 



(6.11) 

(6.12) 



lim 

R w — >oo 



5 - 4tt log 



where the boundary terms are designed to impose (6.8). Note that c m is the lower left 
component of the L m matrix defined in (5.6). Recall that dT>d = U m (O m U O m ) where 
O m = L m {O m ). Also note that the addition of the boundary term at \w\ = R w — > oo is 
to deal with an IR divergence due to the w plane being infinite. Typically one picks the 
generators L m using the freedom to conjugate by PSL(2, C) such that one of the fixed points 
is at oo. Then this IR divergence is absent since w = oo does not appear in the fundamental 
domain. This will not be convenient for us and we choose instead to directly deal with the 
IR divergence 

it was shown that the on shell value of the action (6.11) for (f)$ gives the regu- 



In 



larized action of the bulk gravity solution (up to some minor additions, see (6.13) below.) 
The ZT action captures the conformal anomaly and depends on the choice of metric in a 



fixed conformal class. The appropriate metric here is (6.6) and by the usual dictionary of 



AdS/CFT the action Szt gives a contribution to the partition function Zj^{ds 2 ) of the CFT 
defined on this metric: 



gr 



(n - 2) log A 



(6.13) 



where 7 labels the particular gravitational saddle. We have also included a UV cutoff factor 
oc log A which cannot be removed in the limit A — > due to the conformal anomaly. It is 
proportional to the Euler character of which is — 2(n — 2). 



Putting everything together the partition function we seek (1.2) is the sum over the 
different saddles 7 (including the ones we do not construct): 



Z M (ds 2 ) = J2 e M-SJr + Sl + O(c )) 



(6.14) 



Rather than work directly with Szt and Sl we would like to compute their on-shell 
variation with respect to the Zj. We can use several results in the literature. These results 
can be understood as essentially arising from conformal ward identities. 

For Sl the variation was given originally by Polyakov using the Liouville theory path 
integral. A proof using just the classical Liouville action was given in [33ti36] . The varia- 
tion gives the Fuchsian accessory parameters defined in terms of the Fuchsian monodromy 
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problem: 



dSr 



cn 



a* 6 »> ^ 

For Szt the results in [2H] and [SB] can be applied. Here a mathematically rigorous 
procedure for varying the action with respect to the moduli was used and goes under the 
name of quasiconformal transformations. From the results in these papers we can derive: 



dS 



ZT 



cn 
~6~ 



(Pi ~ Pi) 



(6.16) 



the difference in Fuchsian and Schottky accessory parameters. Details are given in Ap- 
pendix [E} The reason the difference in accessory parameters appears is that the stress tensor 
of Szt which appear when varying the action takes the form: 



T ww = d^ s + -(d^s) 



w 



'(z)- 2 (T zz -T? z ) 



w (z) 




(6.17) 



This stress tensor is the one associated to the ground state of the CFT living on ds 2 



Adding (6.15) to (6.16) explains equation (3.4) given in the prescription of Section 111 



which was the main goal of this current subsection. Note that pf cancels between (6.15) and 



(6.16) so the final result does not depend on the Fuchsian uniformization. This is expected 
since we used the uniform metric ds 2 only as an intermediary. 



B. Regularization Surface 

The goal of the next two subsections is to find an expression for S gr without resorting to 
taking derivatives thereof. To do this we go through the derivation in j27j using a slightly 
different regularization procedure. 

Following [27] we need to pick a regularization surface in order to define the bulk action. 
This surface should be consistent with the L m identifications. Also the desired boundary 
metric should be induced on this cutoff surface . We use the field <ps to define our cutoff 
surface: 

^Ae 0s/2 ds 2 \ A « A- 2 e~* s dwdw = A~ 2 ds 2 (6.18) 
with UV cutoff A 0. Note we could not simply cutoff at fixed £ since under the PSL{2, C) 



isometries (4.2) the coordinate £ changes. For small £ this transformation is consistent with 



7 In the literature on Liouville theory one considers a slightly different form of the Liouville action Sl from 
(6.2 1. Firstly the action is defined on the z plane with reference metric ds 2 . We can get to this form by 



integrating by parts on (6.2) and we go through this in Appendix [C| Secondly the z-plane is not multi 



sheeted, rather the points z% are conical deficit singularities for the uniform metric ds 2 . One can think 



of these as arising from a quotient of our surface M. by Z„. This explains the extra factor of n in (6.15) 
compared to for example |36j . 



30 



the transformation of e^ s ^ 2 . Hence this choice. Using this surface which we call T>a we define 
the regularized action in the usual way: 

167iG N S gr = I d 3 x s /g(R-2) + 2 I d 2 xVh(l — K) (6.19) 
JQa Jv a 

where Qa is the regularized portion of the fundamental domain for the quotient AdS' 3 /Y> 
with boundary T>\. K is the trace of the extrinsic curvature and h is the induced metric. 
Note that depending on the choice of fundamental domain for the quotient Qa the boundary 
will be conformally equivalent to either T> s or T>d as A — > 0. The choice of T>d was worked 
out in [27] so here we pick V s . 

We can be a little more precise and pick coordinates {A, u, u} in order to write the original 
AdSz metric in the Fefferman- Graham expansion with induced metric ds 2 . Following [37] 
write 

4Ae^ s/2 2A 2 e^ s 
e = 4 + A^|W W = U - d ^ S A + A 2 e^ s? ( 6 - 2 °) 

where now (ps is considered a function of 4>s(u,u) which anyway approaches (w,w) at the 
boundary. The bulk metric is: 

2 



, 2 dA 2 1 
ds = A2- + A2" 



1 + A 2 \ 1 



1 - t 2 / 2 



dt + -A 2 (1 - \t\ 2 )T tt dt 



(6.21) 



where we have written the answer in terms of the Poincare disk coordinate t. 

The answer has a particularly simple form since the FG expansion terminates in three 
bulk dimensions [381 EH] ■ The stress tensor T u of the field theory living on ds 2 and in the 
state defined by the saddle at hand appears as a sub leading term in the FG expansion [40J . 
In these coordinates we can write the stress tensor as 8 : 

Tu = (t'(z))- 2 (K z} - {t, z}) = (t'(z))~ 2 (T zz - TfJ (6.22) 



where t, w were given in terms of solutions to the appropriate Fuchsian (6.4) or Schottky 



(4.6) odes . We use the notation T to denote the stress tensor for the vacuum state of the 
theory living on the uniform metric ds 2 . While T is reserved for the stress tensor of the 
theory defined on the singular metric ds 2 . 

Fefferman- Graham coordinates typically develop a coordinate singularity away from the 
boundary. This is indeed the case, the metric becomes degenerate when: 

A 2 = (6.23) 

K-W)\T tt \-Arl 

Actually we could attempt to bypass altogether the Fuchsian uniformization - and never 
even mention ds 2 or the Poincare disk coordinates t. In this case we should pick our regu- 
larization surface such that the induced metric is directly ds 2 which is anyhow the desired 



The normalized stress tensor of the state is related to this by a factor Tcft — 
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metric. The fact that the partition function depends on this regularization surface in the 
limit A — > is a manifestation of the Weyl anomaly in holography [41]. This would leave 
the introduction of the Liouville action Sl unnecessary. We can do this by setting 

= <p s - <p F = lnw'(z) + lniD'(z) (6.24) 

and using to define a new regularization surface. Note that is a locally harmonic function 
on the z plane. Since <pp is single valued on Ai the new field transforms in the same way as 
0s around the non-trivial B-cycles (6.8). 



We can pick Fefferman-Graham coordinates with respect to the Liouville field by re- 



placing 05 — » in (6.20). The bulk metric is then: 

dA 2 1 



ds' 



A 2 A 2 



A 



2 



dz + —T zz dz 



(6.25) 



From this we see that T zz is the stress tensor of the theory living on ds 2 . This is stress tensor 
that appears in the original ode. 

Compared to the Fefferman-Graham coordinates for the metric ds 2 those for ds 2 are rather 
singular. Here we find a breakdown of the FG coordinates arbitrarily close to the points on 
the boundary (z — » z^, A — > 0). This breakdown was discussed and confronted in [42] in a 
similar computation of EREs. They break down at A 2 = 2/\T zz \ ~ \z — Zi\ 2 . This is because 
in addition to UV regulating the theory using the cutoff surface A = constant we need to 
regulate the divergences associated with the conical singularities in ds 2 . We achieved this 
previously by using the singular Liouville field 0^ to transform to the non-singular metric 
ds 2 . Then the Liouville action for (pp contained the divergences associated to these conical 
singularities. So the field has to take into account both the divergences associated to the 
conical singularities as well as the L m identifications. We found it convenient to deal with 
these issues separately by splitting this into two steps. 



C. Action from the Symmetric Domain 



We are now ready to calculate (6.19). Firstly let us compute the bulk integral using 



the (£,w,w) coordinates. We use the FG coordinate (6.20) to define the regulating surface 
at A = const. The fundamental domain was depicted in Figure 11 consisting of removing 
hemispheres from AdS^. Define V the volume of this domain: 



d 3 x^/g(R - 2) = 4V = 4 / d 2 w 



d£ 



d w- 



[w.w) 



(6.26) 



Away from the hemispheres the radial integral is cutoff at £~^ n = A~ 2 e~^ s ( w ' w * > — \d w (j)s\ 2 /2 + 
. . . where we should emphasize that we are working with w,w coordinates at the boundary 
and not u,u. We define V m , V m the volumes of the chunks of AdS^ below the hemispheres 
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segments (below in the sense of Figure 11 ). There are 2n of these but by the replica symmetry 
they are all the same: 

V = \ I d 2 w (2A- 2 e-^ - \d w <p s \ 2 ) + 2nV 1 (6.27) 
4 Jv s 

Where V± corresponds to the "first" hemisphere segment - on the boundary it becomes the 



m = 1 segment of (5.4). The volume is 



where we do the integral using cylindrical coordinates about the center of the hemisphere. 
We have given the radius of the hemispheres p in terms of the Schottky parameter xs defined 
in Section [Vj 



Examining the geometry of the hemispheres shown in Figure [TT] we see that we get two 
terms, one from where the r integral is cutoff by the intersection of the hemisphere with the 
regularization surface and the other from the remaining triangular shaped region. That is 
where the radial integral is cutoff by the intersection with another hemisphere. These two 
terms are: 9 

Vl = ~\i de(<f>s + 2HA/p)) + N 1 , N l = - 1 - jT " d6 In (l - ° ) ( 6 - 30 ) 
The first term is an integral on the AdS% boundary along the segment in the w-plane which 



can be described as (see (5.4)): 



where Wi is the center of the circular arc in the complex plane. Note that along this arc 0g 
is identified under L\ with the 0s at the next arc moving in an anti-clockwise direction on 



Figure 11 We can write this identification (6.8) simply as 



4>s ( wi + r 1 exp(— id + z27r/n)) = <ps (^1 + t exp(i9)) — 41n(r/p) (6.32) 
where W\ is the center of this adjacent arc {w\ = e l27T ^ n Wi). 

9 The expression for Ni is only a function of n and does not analytical continue well to n < 2. We guess 
an expression that has a better continuation in n 

•j / .min(7r/n,7r-7r/n) / 2 ( / \ 

n—^-vl ■«""( 1 -w 1 ) (6 ' 29) 

For n < 2 this expression subtracts the volume of a triangular shaped region, since now the other term in 
V\ over count the volumes of the hemisphere segments. This is a guess since the bulk solution does not 
make any sense for n < 2. This guess seems to yield the correct answer. 
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An issue we have ignored so far is related to the IR divergence associated with working 
in the Poincare patch. To fix this we momentarily move to global coordinates where the 
metric is: 

,2 d^ 2 (R w £ \ 2 2 dwdw 

ds = t + It " r-J R "WTWf (6 ' 33) 

and where the radial coordinate ranges over < £ < R w . The limit R w — )■ oo returns us to 
Poincare coordinates. However before we take this limit we get an extra log contribution to 
to the bulk Einstein action: 

Rw ( 

(Ri + H 2 ) 2 \2% 



V = I d 2 w /m ( + 2 log(UjR w ) +...) (6.34) 



The log term encodes the coupling of 0s to the curvature of the w sphere, which we have 
hidden at \w\ — > oo by working on the plane. We must keep this term which in the limit 
R w — > oo gives us the addition: 

V-+V + - [ d9(f>s - 2ir\n(R w /A) (6.35) 

2 J\w\=R w 



The extrinsic curvature part of the gravitational action ( |6.19 ) is most conveniently eval- 



uated in FG coordinates (A, u, u) - which can be related to (£, w, w) coordinates close to the 
boundary (note if we were using and not <ps as our Liouville field there would be some 
extra complications to deal with here.) That is: 



2/^1^.-/^1^-^) 



(6.36) 



Adding everything together the quadratic UV divergence associated to A vanishes leaving: 

967T - f f 

S gr = - d 2 w ((d(p s ) 2 + Ad 2 <p s ) + ZZnVi + 8 / d6<p s - 32tt ln(i^/A) (6.37) 

c Jv a J\w\=R w 



where one should remove \w\ > R w to define the symmetric domain T> s . The d 2 <ps term 
evaluates to something proportional to the Euler character of M.: 

J d 2 wd 2 (ps = -8n(n - 2) (6.38) 



Combining the gravitational action with the Liouville action (6.2) as in (6.14) we find 
after integrating by parts: 

— S gr = - I d 2 w{d<p) 2 -8n [ dd(<p - 2 In p) + 8 I d6<p - 32tt In R w 

C JV S JU! J\ W \=R W 

+327i(n - 2) In A + 327r(n - 2) + 32nN 1 (6.39) 

where everything is written in terms of = <ps — <Pf- Recall that this is a harmonic field 
that can be defined solely in terms of the Schottky uniformization coordinates w(z) see 
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(6.24). The new domain T> s is a regularized version of T> s denned by cutting out various 
holes where (f) diverges. The justification for this cutting procedure is the same as for the 
Liouville action that we went through in Appendix O Some details have been swept under 



the rug in arriving at (6.39). For example we need to disentangle the conical singularity 
at z 4 from the curvature singularity of the infinite w plane at w = oo. Recall that for our 
choice (5.2) these points were the same w{z^) = oo. The quickest way to deal with this is 
to deform the point w(Zi) — > u>4 1 such that w(z 00 ) = oo for z^ « z^ on only one of the 



replicas. Taking w 4 — > oo of this procedure defines the action (6.39). 



We can now evaluate (6.39) since 4> is harmonic. After some work (the details of which 



are given in Appendix |DJ) one finds an answer which can be succinctly written in terms of ip- 
the particular solution to the ode appearing in the denominator of the Schottky coordinate 
w = \ip+/?p- (see the discussion around (5.2).) We send z\ = 0, z<i = x, z^ = 1 and z^ — » oo. 
By defining the Mutual Information (for this particular saddle 7 = 0) as in (3.26) we can 
take the limit z^ — > 00 without the associated IR divergence. We find: 



?? 



12vr(n- 1) 
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nNi 



+ 



37r(n - I) 

n + 1) ln(n) 
6(n- 1) 



(6.40) 



where ip- ~ z 1 ^ 2 ^ 1 ^ 2 '^ at the origin of the z plane. The integral is along the real axis with 
■0_ = ip-{z + if]) evaluated just above the real axis and we have defined ip- evaluated just 
below the real axis as ip- = ip-{z — irj). Note that ip- is related to ip- by a monodromy 
around the loop on Ai which connects the top of the real line segment z G [x, 1] + irj with 
the bottom z £ [x, 1] — irj. See the left panel in Figure [7j This is the monodromy loop that 
defined the matrix L\ and we can write ip_ = Ciip + + diip_ . Recall that we are studying the 
saddle associated to the monodromy conditions in Tp. A similar expression to (6.40) exists 
for T a . Finally p and /ij are also extracted from tp_ simply as: 



An 



-1 



lim A 



[M 4 



— lim -^k-z 

Z^rOO %h_ 



1-1/n 



(6.41) 



where W is the Wronskian of two solutions to the ode. Note that when we plug these 



constants into (6.40) the factor A drops out so we can effectively ignore it. 



We have confirmed numerically that the expression (6.40) gives the same answer as the 



one obtained by integrating the accessory parameter. We think formula like (6.40) should 



generalize for other non-replica symmetric saddles which would be useful in checking the 
assumption that replica symmetry remains unbroken. The details of course will be slightly 
different and we leave this to future work. 
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VII. DISCUSSION 



To summarize we have calculated some contributions to EREs in holographic CFTs by 
constructing higher genus gravitational handlebody solutions. We only found a subset of all 
possible classical solutions. These are the solutions which were highly symmetric - respecting 
the symmetries of the boundary surface. We found that the bulk actions of this set of 
solutions continued nicely under n — > 1 to the Ryu-Takayanagi formula involving lengths of 
bulk geodesies. We did this for arbitrary numbers of intervals, however the bulk solutions 
were only described in detail for N = 2. 

If we could show that the symmetric solutions dominate in the sum over saddles at large 
central charge for all < x < 1 then we would have found the Mutual Renyi Information 
for all 1 + 1 CFTs with an Einstein gravity dual description and we would have proven the 
RT formula in this case. We have not managed to come up with a proof necessary for this 
purpose and in fact after some thought we are not sure it is true. 

More conservatively for two intervals one should be able to show that our prescription for 
computing the EREs is correct for any n in a perturbative expansion about x = 0. We can 
argue pictorially that the saddle which we called T a dominates over all other gravitational 
saddles in the limit x — > 0. In particular the minimal length of a curve living on Ai (in 
terms of either metric ds 2 or ds 2 ) which is homologous to the cycle C a (see Figure[4| becomes 
parametrically small compared to the minimal length curve homologous to Cp. This is true 
on all replicas. The bulk solution which has the least volume and hence least action will be 
the one where all the short cycles C a are contractable. If any of the other longer Cp cycles 
were contractible then we will clearly get a larger volume subdominant solution. 

Once we have shown T a is the dominant saddle as x — > then at infinite central charge 
the other saddles cannot be seen to all orders in an expansion about x = 0. This leaves 
open the possibility that another non-symmetric solution becomes dominant at some finite 
value of < x < 1/2. The danger region is x ~ 1/2 since then the minimal length curves 
of C a and Cp approach the same length. Similar arguments can be made about the point 
x = 1 as well as for more than two intervals. 

The paper [22] found the same monodromy prescription that we gave in a completely 
different manner. They studied semi-classical (in the sense of large central charge) Virasoro 
conformal blocks for low dimension operators which can be computed in terms of this mon- 
odromy condition [19] . The ERE thought of as a 4-point function of twist operators was 
then argued to receive its dominant contribution from the conformal block for exchange of 
the unit operator (including the stress tensor and its descendants.) These conformal blocks 
can be identified with the bulk solutions we constructed and since they contain the unit 
operator in the s-channel (t-channel) exchange as x — > (x — > 1) must give the dominant 
contribution to the 4-point function about x = (x = 1). The question of intervening 
saddles in the danger region was also unresolved in [22J and would involve the conformal 
block of some other heavy operator (with conformal dimension 0(c)) potentially becoming 
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dominant. 

We now review some material that may help for trying to construct these missing saddles 
in this danger region. 

A. The missing saddles 

Summing over saddles in AdS 3 holography is a well studied subject at genus one (see for 
example [32J S3]) but is less explored at higher genus (however see [IHl 06].) The basic idea 
can be understood in terms of the moduli space of Riemann surfaces M n . At genus (n—1) this 
is a 3 (n— 2) complex dimensional space. The covering space of M n is called Teichmiiller space 
T n and this space distinguishes Riemann surfaces related by large coordinate transformations 
- sometimes referred to as Modular transformations or elements in the Mapping Class Group 
(MCG). M n is then just the quotient of T n by the MCG. If one studies a modular invariant 
CFT on a Riemann surface then the partition function is a function on T n which is consistent 
with the action of the MCG - thus it is also well defined on M n . 

The handlebody solutions to AdS% gravity however are not invariant under the action of 
the MCG. This is because we had to choose a set of (n— 1) A-cycles which were contractable 
in the bulk and these cycles change under the MCG action. This means that the gravitational 
action thought of as a function on T n becomes multivalued on the quotient M n . To find 
a modular invariant partition function we need to sum over all images of the MCG 10 
Schematically, 

Z M = e W (-S^ + O(c )) (7.1) 

geMCG 

where we continue to use a to denote the replica symmetric solution corresponding to the 
monodromy conditions specified by T a . The action g on this solution which we denote 
g(a) then scans through all relevant handlebody solutions. After the sum Z M is modular 
invariant and well defined on M n . 

For example there is one element in the MCG which sends r Q — > Tp . The phase transition 
between these two saddles at x = 1/2 is a fixed point of this action. The other solutions 
g(a) have not yet been constructed. At finite central charge c we would need all of them to 
find a consistent partition function, however at large c we can ignore all but the dominant 
ones. One then just needs to find which is the dominant as a function of x. 

We have found S gr along a one dimensional slice in T n which is the special slice dis- 
tinguished by the replica symmetry. We just need to construct the other one dimensional 
slices of T n related by elements of the MCG to the replica symmetric slice. Since the replica 
symmetry acts non-trivially on these solutions there will be more than one slice of T n space 
(related by the replica symmetry) with the same action. If one of these is dominant then 



not all!, a subgroup leaves the bulk solution untouched, but we do not dwell here on such details. 
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it will clearly give rise to the phenomena that might be called replica symmetry breaking. 
Actually we require infinite central charge otherwise the usual arguments about the lack 
of symmetry breaking in a finite system apply. It would be interesting to understand this 
phenomenon (if it were to occur) in more detail. Some questions that come to mind: What 
is the order parameter? What does this imply about the spectrum of the reduced density 
matrix? What are the implications for the product orbifold theory and EREs written as 
twist operator correlation functions? 11 

Unfortunately finding the other solutions is easier said than done. We sketch how one 
might construct these missing saddles using the methods of this paper. Pick a non-symmetric 



set of [n — 1) A-cycles and demand that they are contractible in the bulk. The ode (3.1) 



can be used to construct the bulk solution - however as we alluded to throughout the paper 



we need to turn on the other accessory parameter (4.9). Which ones we have to turn on 
depends on the symmetry breaking pattern. Already we see the problem becomes intractable 
for large n - we have to search in the 3(n — 2) dimensional space of accessory parameters 
in order to satisfy the monodromy conditions. However it is feasible to attempt this for 
small values of n. One then needs to compute Sgi a \ For this we can no longer resort to 
integrating the accessory parameters and instead have to work with an absolute expression 



for S gT . This is where a suitably generalized version of (6.40) will come in handy. 

To further complicate things there are nonhandlebody solutions when n > 2. See [48J for 
a discussion in the context of partition functions of Ad S3 gravity. See also [32] where these 
solutions are constructed using Quasi-Fuchsian Kleinian groups. Note that they exist if and 
only if the boundary surface has some discrete symmetry, and indeed we sit on a point in 
moduli space with lots of symmetry. We have access to their actions since they are simply 
related to Fuchsian uniformization and the Liouville action Sl [IE] . It is believed that these 
solutions cannot be dominant because they would lead to certain pathologies from the dual 
CFT perspective - relating to studying the CFT on more than one disconnected surface. 
But it remains an open question to show this. 12 We are actually in a position to study this 
question however for now we leave this to future work. 

For the case n = 2 one can construct all the saddles. This is the case where A4 is a torus 
which is then the same surface which is used to compute the thermal partition function of 
the CFT on a circle. As discussed in [20J the transition between the T a and Tp saddles at 
x — 1/2 is related to the black hole Hawking-Page transition. There are also an infinite set 
of saddles [43] which turn out not to contribute at large central charge for a purely imaginary 



ii 



We give an answer to the first question in this footnote: the non-symmetric accessory parameters p s (for 



s ^ 1) are good order parameters, and can be extracted from integrals of stress tensor on A4, see (4.9). 
Formally one would need to slightly break the replica symmetry in order to get a non-zero expectation 
value for this order parameter. Then at large central charge one can remove the symmetry breaking term 

and still potentially arrive at a non-zero answer for (p s ), indicating symmetry breaking has occurred. 
12 We thank Alex Maloney for discussion on this. 
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torus modulus r = 1T2. These come from SL(2, Z) modular transformations of the torus 
and lead to the following saddle contributions to the Mutual Informatiom: 

where A,B e Z and are co-prime. The cycle that is contractable in the bulk is AC a + BCp 
where C Qi p were defined in Figure |4j Note that under the action of complex conjugation 
(C a , Cp) —> (—C a , Cp) so these bulk solutions break this symmetry when both A and B are 
not zero ( in which case there are multiple solutions with the same action.) However as is 



clear from (7.2) the dominant solutions are either (A = 1, B = 0) or (B = 1, A = 0). We 
take this as evidence in favor of the absence of replica symmetry breaking. 

Finally note that understanding how the actions of these missing saddles can be continued 
to n non-integer so we can take the limit n — > 1 is also important. This seems rather tricky 
when the replica symmetry is broken since we have to continue the symmetry braking pattern 
to non-integer n whatever that means. 



B. Further work 

Aside from studying the possibility of replica symmetry breaking we see several avenues 
for how to extend this work. Firstly it is important to understand quantum corrections to 
the RT formula. This involves calculating one loop determinants for fluctuations about the 
saddles that we constructed. 

Secondly one could try to generalize our setup in several ways. For example one could try 
to work with different states in the CFT like finite temperature on a circle. Moving to higher 
dimensions would be difficult because we no longer have the power of AdS^ holography. The 
power stems from the fact that gravity has no propagating modes in three dimensions. Of 
course one could still try to work numerically in higher dimensions, and as we have learned it 
may not be necessary to solve the full numerical problem at integer n in order to reproduce 
RT. It might be possible to proceed by simply setting up the problem well enough so that 
in principle the ERE could be computed. Then if the bulk action can be read off without 
reference to the specifics of the bulk solution it may be possible to continue the answer to 
n — > 1 without doing the hard numerical work. 
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Appendix A: The Riemann surface 



Specializing to two intervals N = 2 the Riemann surface is defined by: 

y n = (g ~ ~ gg) (A1) 

(z - z 2 ){z - z 4 ) 

• z and y can be thought of meromorphic functions on the surface Ai 

• z defines a holomorphic map from Ai to CP 1 branched at the points Zi. The degree 
of the covering map is n. And the branching order is: 

B z (Zi)=n-\ (A2) 

so that the Riemann- Hurwitz relation tells us that the genus of M is g = n(—l) + 1 + 
2(n - 1) = n - 1. 

• The surface has a cyclic Z n symmetry (automorphisms of M) defined by cyclic per- 
mutations of the n different sheets. It is generated by (y — > ye t7T ^ n , z —> z). 

• The space of complex structure deformations of such a surface is 3(n — 2) complex 
dimensional (for n > 2), however we are only interested in a one real dimension slice of 
this space. These will be the deformations which preserve the Z n rotations discussed 



above and leave the surface in the form described by (Al). 

A basis for the holomorphic differential forms is: 

dz 

v l = — , t = 1, . . . , n - 1 (A3) 

y\z - z 2 ){z - z A ) 

A linearly independent basis for the quadratic holomorphic differential forms is: 13 

dz 2 

" s+1 = — w w w r s = 0,...,n-2 (A4) 

y s {z - zi)(z - z 3 )(z - z 2 ){z - Zi) 

dz 2 

u s+n ~ 2 = — - 2 s = 2,...,n-2 (A5) 

y s (z - Zi)(z - z 3 )(z - z 2 y 

dz 2 

^ +2n " 5 = — { v Tt Vi s = 2,...,n-2 (A6) 

y s (z - z\){z - z 3 )(z - z A ) 2 

d7 2 

u 3{n ~ 2) = — (A7) 

There are 3(n— 2) of these, and they label the space of complex structure deformations. 



13 



In order to find these one notes that the point defined by z — > oo is smooth and thus u> takes the form 
P(z)dz 2 y~ s where P(z) is meromorphic function with poles at the branch point and net pole order —4. 
There are linear relations that need to be taken into account between these meromorphic functions. Then 
the bounds on s can be found by examining the behavior close to each branch point. 
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• Only one of the quadratic differentials is invariant under the Z n symmetry: 

1 — ^z 2 

(z - zi)(z - z 3 )(z - z 2 )(z - z 4 ) 

As usual we can use an SX(2,R) transformation on z to move the points ^1,^2,^3,^4 to 
0, x, 1, 00 respectively. After scaling y by a constant the surface can be described by: 

y" - f^l < A9 > 

Appendix B: The connection matrices 

In this section we setup the connection matrices which allow us to study the monodromy 
problems numerically. For each point Zi consider the two linearly independent solutions: 



where z is taken slightly to the right of Zi on the real axis. Since we know the ode is real 
along the z axis it is useful to construct the following real connection matrices: 



u 



(z) = R i;i+1 Tu^\z) (B2) 



where T enacts the monodromy of moving z from slightly to the right of Zi + \ to slightly to 
the left of Zi+i 

T -{ -ie*"™) (B3) 

From this it is clear that the remaining matrix -R^+i is real. Further the Wronskian condition 
on two linearly independent solutions tells us that det Ri^+i — — 1. 

To summarize we have defined the following set of 2N real connection problems: 



i+1 



if> « a i+1 (z i+1 - z)5+2 1 n + p i+1 (zi +1 - z)* 2 «, z->zl 



(ai+i A+i) = (at ft) R 



(B4) 

which are amenable to numerical work. From these i?'s and the matrix T we can construct 
the full set of monodromies in a straightforward manner. For example the trivial monodromy 
condition on a path encircling two adjacent points (zk,Zk+i) is: 

1 = T 2 R k>k+1 T 2 R-* k+1 (B5) 
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Appendix C: Regulating the twist operators 



The Liouville action (6.2) needs to be regulated by smoothing out the singularities of 4>f 
close to the singular points Zj. Following [37J we cut out infinitesimal holes of A4 around 
these points \z — Zi\ < e and "insert the unit operator" by picking a different metric within 
these holes: 

ds 2 = dzdz (\z - ^| 2 /e 2 ) 1/n_1 oc dtdie 2 - 2/n \z - z { \ < e (CI) 

and we have matched onto the metric ds 2 = dzdz at the boundary of the hole. Doing this 
allows is to extend the definition of (ftp to the closed surface Ai. Note that for this choice of 
metric <pF is essentially constant in this hole (up to corrections due to the curvature of the 
hyperbolic space.) There are 4 of these holes. 

Additionally we also have to chop out z > R z which contributes an IR divergence to the 
path integral. There are n images of these z — > oo holes in the t plane at t = t™. To do this 
we work with the following metric for large z: 

n zn z — 

ds 2 = R 4 zT1I - oc dtdt, z> R z (\t - CI < #/R,) (C2) 
which again matches to ds 2 = dzdz at the boundary and sets <pp to approximately a constant 



in this region. In this way the action Si in (6.2) is well defined. It is easy to show that the 



contributions to Sf from the n + 4 holes discussed above is zero essentially because 0^ is 
a constant and the volume of the holes is going to zero in the limit e — > 0,R Z — > oo. We 
can then restrict ourselves to performing the integral in the Liouville action over the region 
outside the holes: M = M\{(\z - z { \ < e) U (\z\ > R z )}. 

We now write this integral in terms of the original z coorindates: 

8 L = ^- / d 2 z ((«9<M 2 + 2<M<9 2 <M) 

96tt Jm 3 

where we have integrated by parts and used the fact that the Euler character of the surface 
Ai is 2(2 — n) to introduce the Liouville term oc e~^ F into the action. Up to this constant 
the final result is that of the standard Liouville action for the field 4>f living on the z-plane. 
The boundary term also taking the usual form: 

(—) s hdry = -4 ( i - -) [ de^ + s V / de^ (C4) 

\ c / V nj J\ z - Zl \= t ^J lzl=Rz 

where the last sum is over the different sheets of the branched covering. Also the first 
integral above is from 6 = 0, 2im. We have used the behavior of <pF close to Zi and z = oo: 

1 - - ) In \z - Zi\ 2 + . . . , z ^ Zi\ 4> F ~ 2 In \z\ 2 + . . . , z -> oo (C5) 
n I 
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This action actually defines the problem that we have already solved. That is if we vary the 
action we find the Liouville equations of motion: 



d 1 



-Se 



(C6) 



now subject to the boundary conditions (C5). There is a unique solution to this problem. 
That is the solution which defines a constant negative curvature on the space Ai and which 
is equivalent to the solution found using the ode and Fuchsian monodromy condition. 



Appendix D: More on the symmetric domain action 



We continue the calculation of the bulk gravity action starting from (6.39). This action 



is integrated over the regulated domain D s defined analogously to Ai in (C3). That is 
D s = D s \{(\z - Zi\ < e) U 0| > R z ) U (\w\ > R w )}. 

The field appearing in (6.39) <p is the Weyl factor for the the conformal transformation 
from the w to z planes: 



ds = dzdz 



'^dwdw 



(Dl) 



Thus it will have singularities whenever either of these metrics is singular (unless the two 
metrics have the same singularity and they cancels out.) It is clear that <fi is harmonic 



so we can evaluate the bulk term of (6.39) by integrating by parts. This then reveals the 



singularities located at zi and as well those at z — > oo which we remind the reader are 
actually n points one from each replica - they are located on the w-plane at w;™. We also 
still have the w plane singularity at w — > oo. Altogether we have: 



1 1 1 71 

l S gr = -<P(\w\=R w )--lnR w + -J2<P(\^\=Rz{w = wZ}) 

"X>(l 



n 



m=l 

n 



Z - Zi 



1=1 



24tt 



(Ji - 8iVi) + — — — In A 



(n-2) 



(D2) 



where we have defined the integral: 



h 



dd (0-41np) 



(D3) 



and the notation is for this integral the same as in (6.30). 



In order to make sense of (D2) we have temporarily deformed the space w such that 
z^ no longer maps to w = 00 but instead maps to w 4 3> 1 and also such that there is a 
point Zoo ~ (on one of the replicas) which does map to to w = 00. This will disentangle 
the w-plane curvature which is hidden at w — > 00 from the twist operator at z = z^ which 
previously mapped to w = 00 (see Figure [7]). The behavior of can then be read of from 
the behavior of the Schottky uniformization coordinate w(z): 
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Around the twist operators z — > zc 

w(z) « Wi + fa(z - Zi) l/n ->■ 0^ 2 ( - - 1 J \n\z-Zi\ + 2\n\fa/n\ (D4) 



• Around the w-plane curvature singularity at w — >• oo (recall that ^ is close to z^ and 
only on one of the replicas): 

w{z) « — — — — >■ « 41n |w| — 21n (D5) 

• Around the z-plane curvature singularities at z — > oo 

iu(z) « w™ + ^ ->■ « -41n|z| + 2^1^! (m = l,...n) (D6) 

z 

where we have defined many new constants fa, [i™,^ which can be extracted from w(z) 



numerical. Plugging this into (D2) we find: 



4 „ „ n-1 



c- l S gr = -^-^^ln^/nl-^lnl^l + ^lnl/i-l (D7) 



12 ^ Ir ' 1 6 11 6 

i=l m=0 



(n — 2) n , , n , fn — 2) , (n — l) 2 , 

-8iVi) - -kiRz- -lnA+^ Mne 



3 24tt v ' 3 3 " 3n 

Note that \/J,™\ will actually all be the same by the replica symmetry. We now want to return 
to the setup of interest (after the aforementioned deformation of the w plane) by taking the 
limit w{zi) — > oo and thus z^ — > z±. We also want to eventually take the limit z± — > oo 
(and set Z\ = 0,Z2 = x,z% = 1) since this is most convenient for numerical work. We take 
these limits sequentially: 

• In the coincidence limit close to w ~ W4 ~ 00 one can argue that the analytic map 
has the form: 

z « z 4 + wffa n (wl l - w- 1 )" 1 (D8) 

which reproduces the behavior about w = W4 and w = 00 if we additionally have 
= —nw2 +1 /n2- Expanding (D8) for large we have: 

w fa{z — -24) _1//n , fa = —w\jfa (D9) 

where we should fix fa in this limit. The we find: 

hm i^r 1 ! = |/r 4 | (n+1 V M (D10) 

• Similarly in the limit — > 00 we have: 



lim |(^) 2 >I n_1 | = |Air 1 tf B " VB >- 2/n (DH) 

24— >oo 



{ 1 j2n„-n-l| _ | 77 n-1^2(n-l/n)| /n _2/ 

where we have defined as z — > 00: 

w « faz 1/n (D12) 
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Taking these limit in (D8) gives: 



c^S, 



12 



In 



/' ; 



" (/l _ 8Nl) + (jlJI Hn) 



in 



2) 



247T 



(n 2 - 1) , 



{n-lf 



lne 



n 



lni? 2 



6 

(n ■ 



(D13) 



2) 



In A 



6n 3n 3 3 

To deal with all the regulator factors we simply compute the mutual information. We need 
the ERE for a single interval using the same technique as above, so that the non-universal 
pieces cancel. One finds (the uniformization map can be found in this case analytically, we 
do not go through the details which can be for example found in [4"T]): 



-1 qN=1 



n 



1 In 1 

ln|zi-z 2 | Inn In + - In A (D14) 

3 3 3 3 



where Zi, z<i are the end points of the single interval. We also need the partition function Z\ 



for the theory on a single replica in order to relate Zm to the EE. See (2.5). This can be 



computed as above to find \a.Z\ = c/3(ln(_R 2 /A) + 1). Putting everything together we find: 



12 



In 



^1/^2/^3 



/'I 



+ 



n 



2Aix(n - 1 



-(/ 1 -8JVi) + 



6 



1 \ , . . (n + 1) ln(n) , 
- + 1 In x + 1 - ] ; D15 
n J o(n — 1) 



where we remind the reader this computation was for the r g saddle. After mapping the 



integral l\ in (D3) to the z plane we get the expression quoted in (6.40). 



Appendix E: Quasi-conformal transformations 



A variation of the modular parameter of the surface Ai can be thought of as a qua- 
siconformal transformation. A transformation which deforms infinitesimally the complex 
structure of the manifold Ai. It was shown in [2H] that the ZT action behaves nicely under 
quasiconformal transformations. The results make physical sense since the variation gives 
an integral of the stress tensor associated with the Schottky uniformization: 

1 



5S 



ZT 



24tt 



M 



d z w ( SNIT WW 



yj-L WW 



T 



:(d w <fis) +d w d w (p s 



(El) 



The quasi-conformal variation is defined through a Beltrami differential <5A^.To describe 
this consider varying the complex structure of Ai by defining new holomorphic coordinates 
on the different coordinate patches. For example working on the branched covering: 



z' = z + 5z(z,z) y' = y + 5y(y,y) 



(E2) 



and equivalently for the antiholomorphic coordinates. Demanding that the new transition 



functions (defined by the complex curve (Al)) are holomorphic with respect to the coordi- 
nates (z',y', . . .) defines a Beltrami differential: 



8N% = d- z 5z(z, z) 



(E3) 



45 



which transforms as the placement of the indices suggests. It is this type of differential that 
appears in (El). The integral in (El) is simply over the entire surface Ai without boundary 
(compared to the integrals over the fundamental T>i) since the action of S does not change 
the stress tensor T and 5N is well defined on the whole surface. 
Let us transform this result to the z plane: 



5S 



ZT 



C 

24^ 

c 
24^ 



M 



d 2 z { T zz d z 5z + Tg S dJz 



f zz = {w,z} - {t,z} 

{pi - Pi) 



V / d 2 z (d z 5z(z, z) ^ 1 P 9 + d z 5z(z, z 



[Z - Zi 



(E4) 
(E5) 



where we have replaced the the Schwarzian derivatives with the their expansion in terms of 
accessory parameters (4.9 ). In particular we have set all replica symmetry breaking accessory 
parameters to zero. We have picked Sz to be the same on each replica and in this way the 
quasiconformal transformations leaves us on the moduli space of Riemann surfaces defined 
by flAlb. 



Note that what we have arrived at in (E5) looks like a Ward identity for conformal 
invariance. Integrating we find, 

^2((Pi -pf)Szi + c.c.) 



5S 



ZT 



12 



(E6) 



where 5zi = 5z{z il Zi). We have integrated over all the replicas explaining the factor of n. 
A certain amount of regularity in 5z(z, z) was assumed in order to be able to invert the 
operator d z in (E5) and this can be justified by giving an explicit expression for 5z in terms 
of 5zi 

5z(z, z) = - l -e^ Sz ^A<Pf (E7) 

i 

where 0^ is the Fuchsian Liouville field and this equation is the same on each replica. See 
for the complete discussion. 
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